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The approach to the calculation of quantum dynamical correlation functions is presented in the framework of the 
Mori theory. An unified treatment of classic and quantum dynamics is given in terms of Weyl representation of operators 
and holomorphic variables. The range of validity of an approximate quantum molecular dynamics is discussed. 
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Methods that permit to determine, at least approximately, the thermodynamical properties of quan- 
tum dynamical systems by means of calculations of classical type usually turn out to be very useful. 
For the static case, this goal has been successfully achieved by means of a variational procedure based 
on path integral formulation of statistical mechanics Q . The method leads to the construction of an 
effective potential that accounts completely for the quantum corrections of the quadratic part of the in- 
teraction (the harmonic oscillators) and treats in a self-consistent one loop approximation the quantum 
effects due to nonlinear terms H, [| f|] . The effective potential has proved to be a very powerful tool for 
concrete computations j^, ^j, [t]] : in particular it has allowed for the reduction of quantum Monte Carlo 
to classical Monte Carlo numerical calculations , with evident fruitful consequences || . 

A rigorous treatment of the dynamical problem is, as always, much more difficult. A success in 
such a program, however, is scientifically very remunerative. Indeed, for instance, it opens the way to a 
sound formulation of molecular dynamics for quantum systems. Using the effective potential obtained 
in J^, [| , some naive attempts of extending the treatment to a dynamical setting have been tried |j] . 
However, by admission of the authors themselves, these attempts assume some ad hoc multiplicative 
relation in order to connect the "centroid molecular dynamics" to the quantum correlation function 
ill- 

In this paper we shall present a new more rigorous approach, based on the Mori theory and on 
projection techniques |pj|. At lowest order, but within the framework of a consistent analysis, we 
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recover the heuristic expressions of ref. || showing their limitations and indicating the correct method 
to improve them. 

The key points of our theory are: i) the dynamics of classical and quantum oscillators is the same; 
the only difference is the static amplitude factor determined by the initial conditions, ii) the use of the 
Weyl representation p3[ which translates quantum problems into a classical framework and permits a 
unified approach by which any semiclassic expansion is obtained in natural way fL2f . 

Starting from what we have learned in the theory of the effective potential |], we develop a 
semiclassical calculation in the Weyl representation that describes quantum mechanics in terms of 
the classical canonical variables q and p : expressed in terms ofthe holomorphic variable z = 2~ 1 / 2 (q—ip), 
z* = 2 _1 / 2 (<7 + ip), corresponding to the creation and annihilation operators a) and a. If we consider 
a dynamical system with a Hamiltonian H, the quantum equation of motion of an operator A, in the 
Weyl representation, reads |y| 

^L = ln w sbx^A w =iCA w . (1) 
dt h 2 w 

Here A is the differential operator that gives the Poisson bracket, 

A = i (d z ~3 Z * — *d z * ~3g] . (2) 

and C is called the Liouville operator. A w and TC W are the Weyl symbols of the corresponding quantum 
observables. Explicitly 

0™(z*,z) = e- zz * [^e(-^ + **<-<* z )(z-t\O\z + 0, (3) 
J ni 

where \z) is the coherent state, defined as 



\z) =e\ z * z e z * a] ~ za \Q). (4) 

The canonical equations are obviously recovered when h — > 0. Moreover the statistical average of a 
dynamical variable reads 

(O) = Tr (p O) = J ^ p w (z* , z) O w (z*,z) (5) 

and, again, the classical limit reproduces the phase-space average. 

For sake of simplicity, we shall only consider a system with one degree of freedom and with a 
Hamiltonian 

H = hajQ z*z + gV(z* 1 z), (6) 
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where g is the quantum coupling. In the following, when not necessary, we put h = l. The equation of 
motion for z and z* can be put into the following form along the line of the projection- Zwanzig-Mori 
theory Q: 

dzi ^ = iujz(t) + [ dr M(r)z(t -t) + f w (t) . (7) 



dt 

A similar equation can be written for z*(t), that, as a quantum observable, must be treated indepen- 
dently of z. 

This equation can be solved once the initial condition z(0) is assigned. In the classical context one 
has to fix the values q(0) and p(0), while in the quantum case the initial state of the Hamiltonian must 
be specified(^). In the following we shall use the shorthand notation z(0) = z and the analogous for 
other dynamical variables. 

The term f(t) describes the effect of the fluctuating forces, caused by the non linear interaction 
gV(z*,z) and has the following expression: 

/"(t) =e it ( 1 - p ^ C (l-P)£z, (8) 

where P is the projection operator onto the direction of z. The explicit expression of P, as well as the 
memory kernel M{t) are conveniently given in terms of a scalar product introduced by Mori [TT| . The 
latter plays an important role in the theory and, for quantities with vanishing thermal averages, reads: 

(A\B{t))= f dX((A w T(0)B w (t + iX)). (9) 
Jo 

It turns out that: 

P A B(t) = (A\B)- 1 (A\B(t)). (10) 

One realizes that this is just the Kubo function and in the classical limit reduces to the correlation 
function up to a multiplicative factor. 

The quantities of physical interest are the Fourier transforms C(u>) and R(w) of the correlation 
function, (a'a(i))- probed by scattering experiments - and Mori product, (a\a(t)) respectively. They 
are connected by the "detailed balance" principle: 

1 - e 

The frequency u> is given by: 

iuj = {zlz)- 1 (z\PCz) . (12) 
while the memory kernel M(t) turns out to be 

M(t) = (z\z)- l (f w (t)\f W ), (13) 
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with f° = f w (0). Referring to eq.(Q), M(r) can be written in the form 

M(t) =S 1 S 1 , S x = (fflftTHfiWf?) (14) 

where /q" — z, loq — u, Pq — P, Cm) = C and 

6 1 = (z\z)-\f?\fn- (15) 

It will be necessary to introduce some further definitions in order to describe the recursive framework 
of the dynamical theory. The new quantities at different orders will be denoted by different indices. We 
want, however, to stress that the factor Si is defined only in terms of static quantities. 
The equation of motion for 

E (t) = (z\z)-\z\z(t)). (16) 

becomes: 

^ = iwEo(t) + Si j drEi(T)E (t-T). (17) 

Denoting by Pj a family of projection operators in the direction of the vectors fj, from ( p"7| ) we 
easily generate a hierarchy of equations with reduced Liouvillians: 

Cj = (1 - Pj - Pj-i - -Pi- P )Co . (18) 

Finally, the continued- fraction expansion for the Laplace transform So(z) = R(z)/R(t = 0) is obtained 

i ( f w \ r ( -f w I f w \ 

aj(Z) ~z + iu j + 5 j+1 Z j+1 {z)> Uj ~ {ff\ff) ' d]+1 ~ {ff\ff) ■ U9j 
The quantities u>j and <$j+i are related to the first 2(j + I) moments of the time series expansion 

of So(t). While (/j"|/j") with j ^ can be expressed in terms of static correlation, the normalization 
quantity (fo\fo) requires the direct evaluation of the Mori product, giving the static susceptibility. 

We observe that the calculation of S(z) - and therefore of C(u>) - can be approached from the 
knowledge of the static quantities Sj up to a sufficiently large number j = J , supported by some 
insight into the long time behaviour of the continued fraction termination . The effective 

potential theory permits its approximate determination JtJ and recently the averages of the imaginary 
time ordered products of operators that involve both momentum and displacement have been evaluated 
& 

Let us consider the parameter / = f3Tiuj/2, introduced by Feynman, which determines the quantum 
character of the system through the Langevin function: L(f) = (coth f — 1/f) , From the results of 
|]l6f , we derive explicit expressions for the average of an operator O 

I j dz*dz {{OW{z . z)))e -(rHeff(z*,z) j (20) 



4 



and for the Mori product, at lowest order in /: 



(A\B) ~ ±fij ^^{(A w (z*,z)})((B w (z*,z)})e-P H eff( z *> z ) , (21) 

where TC e ff is the effective Hamiltonian, ((O w (z* , z))) is the two-dimensional gaussian average of the 
Weyl symbol of the operator 0(p, 0] and 

r d£dz_ e -f3H eJ1 {z*,z) (22) 

is the partition function. 

For instance, by this expression, ( ^o|) the static average of a quantity dependent on the coor- 
dinates only, turns out to be 



W-jfijj?f«MQ))'-f nr «f«>, (23, 

where 

((A(0» = pcwr)" 1 / ^A(e + ry) e ^ 2 /(2«) (24) 



denotes the gaussian average of the variable A(x) with moment a — %/ (2mu>) L(f) . 

For A = q 2 the relation ( f23[) gives (g 2 ) = a + h/(2mfu>), where the first term a is a consequence of 
cq.(^4|), while the second term is due to the "classical" average with the effective potential. This means 
that the gaussian average described in eq. ( pi| ) is essential to yield the correct quantum occupation 
number (era) from the classical one. This is a crucial point for a good approximation of quantum 
dynamics. 

A perturbative approach of the Mori theory can be performed when the Liouvillian C can be split 
in two terms C d and C nd satisfying the relations: 

PvCPo =C d , (1 - P Q )CP Q = C nd (25) 

i.e C d cannot connect diagonal with non-diagonal terms and viceversa, while C nd makes always this 
connection. In an exact splitting of the operator C we should extract all diagonal contributions from 
the non linear terms. This can be done, by a one loop approximation, which causes renormalization 
of frequency and coupling. This renormalization is caused both by quantum and thermal gaussian 
fluctuations. 

Here, we want to propose a different perturbative scheme, in the spirit of the effective potential 
method. We account exactly for the classical part as well as for the fully quantum contribution of the 
harmonic oscillators. The purely quantum effects of the nonlinear terms will then be calculated in a self 
consistent way at the level of one loop approximation. In particular, when operating on the nonlinear 
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part gTi-nii a t lowest order in the non linear quantum coupling we shall be allowed for the following 
substitution: 

{2/%)gHZ sin(M/2) ~ gH^A. (26) 
Accordingly, we split the Liouvillian operator into the sum C — C h + gC an , where 

£' l = 2ft(z*z+i)sin^, gC an =gW(z*z)k (27) 

The parameters ft and W will be determined by physical requirements. It is apparent from Eq. (p7j) 
that the quadratic part is treated in fully quantum way, giving the quantum evolution of the "best" 
harmonic oscillator, while the non linearity causes a "classical-like" dynamics. As a consequence of 
the splitting £ = £ h + £ an , at the lowest nonvanishing order in the coupling constant, we have these 
important consequences. The quantum effects on dynamics refers to the imaginary part of the path 
(0, i(3). We take into account that the parameters of the effective potential are calculated at the average 
point of this path along which it is affected by an harmonic potential. As a consequence, fi and W are 
exactly the same as defined by the effective potential for the static case. 

The validity of the use of the effective potential is related to the quantum coupling g. At the first 
order, corresponding to a range of temperatures, f3 < l/g, the definition (JsJ) of fi(iX) becomes: 

/f(tA) ~ e ~X(l - P )£ h fw = e -X£ h (i_p ) g £a» 2 . (28) 



With the approximate expression (|28|) of the fluctuating force, the integration on A in the relaxation 
function can be easily performed. The result of this integration, is cancelled out by the denominator 
in the expressions of 3j,, so that only correlation functions done by the rules given by eqs. (|20|), 
(|23| ) are involved and the parameter Si is calculated with the same receipt. The procedure can be 
iterated at successive levels, using again the approximate Liouvillian of eq.(p7j). Therefore a continued 
fraction expansion is obtained where each coefficient <5's is calculated within the afore-mentioned scheme. 
However it is worthwhile to notice that the nonlinear quantum effects on <5's are evaluated only at one 
loop, so that their inaccuracy is increasing with the order of S n . 

Furthermore, we can do a rather crude approximation for the memory function M(t). Let us 
neglect the Gaussian quantum spread ( p4[ ) of the product of the fluctuating force, under the condition 
ifi(t)fi) w ~ fi(t)fi- This corresponds to neglecting the quantum effects of non commutability of the 
operators in the expression of the memory function and using the classical expressions for occupation 
numbers. 

At this level the equation of motion for z(t) becomes 

d ^±=i£lz(t) + Si f drSi(r)z(t-T)+/r(f), (29) 
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where Cl is the frequency Q renormalized by the "diagonal" part of the classical non-linearity JTl|| , 
containing only correlation functions, (see Eq.[l^) . At the same level, M(r) = (5i2i, is given simply by 
a classical- like correlation function. Therefore eq.(p9f) represents a "classical" evolution of z(t) where 
the frequency f2 and the potential W must be determined in such a way to reproduce the same static 
average values of the observables. At this approximation level, therefore, they must coincide with the 
quantum-renormalized frequency and with the effective potential parameter of the static case 0|. At 
the same level the Mori product 3o(t) is nothing but thedynamical correlation function calculated by 
the classical-like effective potential. Finally, through the detailed balance principle, we can recover the 
spectral shape as guessed in ref. § . 

Of course this approach is valid for g — » 0, i.e. classical case or non interacting harmonic oscillators. 
As the coupling g increases, its validity is limited to short times. The range is dependent, firstly on the 
smallness of the parameter /, because we have neglected the gaussian quantum spread of the fluctuating 
forces, secondly on the inverse of the quantum coupling g. Indeed, only the normalized second moment 
Si is rather well reproduced. All effects related to line shape, like finite lifetimes of the elementary 
excitations are considered at classical level so that there is a shift of the peak of the lineshape, without 
any substantial modification of the linewidth, Jl7[ | - 

A further improvement, proposed by can be suggested, from an heuristic point of view, by the 
lowest order expression for the Mori product in the memory function, |]l6| shown in eq.(pl|). At this 
order the Mori product can be reduced to the average, done by the effective potential, of the product of 
the gaussian quantum averages of the fluctuating forces. This observation could justified the insertion 
of the gaussian quantum average of the forces in the "centroid molecular dynamics" ((^, 10 ). 



Our derivation shows the limitations of a "quantum molecular dynamics" done by means of the 
simple use of the effective potential in the classical equations of motion. The simple corrections of 
the occupation numbers, recently proposed ]9|, [ic|| , cannot satisfy the request of the calculation of the 
dynamical correlations for quantum-interacting (strong correlated) systems. The necessity to account 
for quantum effects in the life-time of the excitation would require anharmonic corrections to effective 
potential on the same kind of the ones proposed in heuristic way j|, |l8). Nevertheless the excellent 
results for the thermodynamic quantities []l8[| this correction has not been completely justified up to now 
nor its extention to the static quantities appears to be straightforward. The quantum computations 
of dynamical correlations is still an open problem which has to be faced by different approaches. The 
moment calculations and the consequent expansions fj], ||] although not exhaustive, are still competitive. 
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